On-Chip Microwave Fock States and Quantum Homodyne Measurements 
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We propose a method to couple metastable flux-based qubits to superconductive resonators based 
on a quantum- optical Raman excitation scheme that allows for the deterministic generation of sta- 
tionary and propagating microwave Fock states and other weak quantum fields. Moreover, we intro- 
duce a suitable microwave quantum homodyne technique, with no optical counterpart, that enables 
the measurement of relevant field observables, even in the presence of noisy amplification devices. 
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A two-level atom coupled to a single mode of a quan- 
tized electromagnetic field is arguably the most fun- 
damental system to exhibit matter-radiation interplay. 
Their interaction, described by the Jaynes-Cummings 
(JC) model, arises naturally in the realm of cavity quan- 
tum electrodynamics (CQED) in the microwave [l| and 
optical domains 2]. There, a variety of nonclassical 
states (e.g., squeezed, Schrodinger cat, and Fock states) 
and other remarkable phenomena and applications (e.g., 
entanglement and elements of quantum logic p() have 
been proposed and realized. Other physical systems, like 
trapped ions 4], can reproduce the JC dynamics and, 
consequently, exploit this analogy for similar purposes. 
The intracavity field in CQED and the motional field 
of the trapped ions are typically detected through a suit- 
able transfer of information to measurable atomic degrees 
of freedom [l|,|4( . In the case of propagating fields, quan- 
tum homodyne detection p , a method based on efficient 
photodetectors, leads to quantum state tomography |g. 

Given its relevance for quantum communication, 
single-photon sources have been pursued in the optical 
domain 0, 0. Recently, several CQED-related experi- 
ments have been performed in tunable, solid-state sys- 
tems. Quantum dots in photonic band-gap structures 
have been used as single photon sources p and super- 
conducting qubits have been coupled to on-chip cavi- 
ties UjJ- In addition, microwave squeezing with 
Josephson parametric amplifiers [l2| and aspects of the 
quantum-statistical nature of GHz photons in mesoscopic 
conductors have been demonstrated 0. 

In this Letter, we propose the implementation of a de- 
terministic source of microwave Fock states, or linear 
superpositions of them, at the output of a supercon- 
ducting resonator containing a flux qubit. A Raman- 
like scheme |l4| determines the coupling between the 
cavity and the qubit that consists of two metastable 
states of a three-level system in a A-type configura- 
tion 0, EH Ell EH- Furthermore, we show that rele- 



vant observables of these weak quantum signals can be 
characterized by a microwave quantum homodyne mea- 
surement (MQHM) scheme. This is a speciall y e ngi- 
neered technique based on a hybrid ring junction |l£j,[2£| 
acting as an on-chip microwave beam splitter (MBS). 
Given that direct microwave photodetectors are not 
yet available, the measurement device consists of linear 
phase-insensitive amplifiers followed by square-law detec- 
tors (and/or mixers) and an oscilloscope (or spectrum 
analyzer). These fundamental tools represent building 
blocks to establish on-chip quantum information trans- 
fer between qubits. 

A prototypical example of a flux-based quantum circuit 
is the radio-frequency (RF) superconducting quantum- 
interference device (SQUID), a superconducting loop in- 
terrupted by a single Josephson tunnel junction. The RF 
SQUID Hamiltonian is 
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where Q is the charge stored on the junction capacitor Cj , 
4> is the total flux threading the loop (with [<£, Q] = ih), 
<£ x is an externally applied quasi- static flux bias, L s is 
the self-inductance of the loop, Ej = I c o^o/2tt is the 
Josephson coupling energy, I c o is the junction critical 
current, and $o = h/2e is the flux quantum. For appro- 
priate design parameters, and close to half-integer values 
of <£ x /<£o, the RF SQUID potential profile becomes a 
relatively shallow double well whose asymmetry can be 
tuned by setting <£ x (see Fig.GJ. In this case, the two low- 
est eigenstates |g) and |e) are metastable and localized in 
the left and right wells respectively, whereas the excited 
state |h) is delocalized with energy above the barrier. As 
will be shown later, these levels are suitable for imple- 
menting a Raman excitation scheme. The energy levels 
can be tuned by statically biasing <£ x during the experi- 
ment and transitions between states are driven by pulsed 
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FIG. 1: (Color online) (a) Potential profile of the RF SQUID 
for two distinct values of <3> x [(I) and (II)] and wave func- 
tions for the first three lowest energy levels (|g), |e), and |h)). 
Parameters are given in Table I. The Raman scheme is indi- 
cated with arrows, (b) RF SQUID energy-band diagram near 
$o/2 plotted vs. 4> x . Zero detuning case, point (I), and large 
detuning case, point (II). (c) Absolute value of the vacuum 
Rabi frequency ghe as a function of <J> X - The coupling reaches 
56 MHz at the anticrossing between levels |e) and |h). 

AC excitations. Because of their relatively large excited- 
state lifetimes, persistent-current (PC) qubits [ll| can 
also be considered for the purposes of this work. 

The segment of superconducting coplanar waveguide 
(CWG) shown in Figs. (a), (b), and (d) is one real- 
ization of an open-circuited transmission line resonator 
capacitively coupled to a CWG transmission line. Such 
a cavity is characterized by eigenenergies with transi- 
tion angular frequencies uok — 2nfk that are much larger 
than the thermal energy at cryogenic temperatures. 
Its Hamiltonian is H c = J2k (^k^k + V^) ■> where a\ 
and a k are the bosonic creation and annihilation oper- 
ators for mode k. Current and voltage, corresponding 
to magnetic and electric fields respectively, are conju- 
gate operators associated with the quantized resonator, 
I c (z,t) = (d/dt) , &(z,t), where z is the spatial coordi- 
nate for the superconducting inner strip and i)(z , t) is 
the normal mode expansion of the cavity field [2l|. The 
vacuum current for the odd, fc Q , and even, fc e , cavity 

modes is I^ ko/e ( z ) = Y.k° =i V^^TC^O cos [k nz/D) + 

J2k=2 V h^ke I (Dl) sin (k e 7Tz/ D) , where k , K e are up - 
per cutoffs for the odd and even modes respectively |2l[, 
D = A/2 is the length of the resonator (A is the full 
wavelength corresponding to (Jk =i), and I is its to- 
tal series inductance per unit length. Hereafter, the 
cavity is chosen to be operated at fco = 1 and is as- 
sumed to have an external quality factor Q x ~ 10 4 at 
fi = ooi/2tt ^ 10 GHz, corresponding to a cavity decay 
rate n c ~ 1 MHz Hence, at a base temperature 

Tb ~ 50 mK, the mean number of thermal photons is 
(^th) = [exp(/kJi//cBTb) — « 10 -4 and the cavity 
mode can be considered to be in the vacuum state |0). 
Embedding the RF SQUID in the CWG cavity allows 



an inductive coupling between any two levels of the RF 
SQUID and a single cavity mode k. The resulting inter- 
action Hamiltonian is i^j = —M CS I C fc J s , where I Cj k and I 8 
are the resonator and the RF SQUID current operators 
respectively and M cs is their mutual inductance. Using 
the explicit forms of the current operators, we find 

H I = - (M cs /L s ) I° C}k (z) ($ - * x ) i [a{(t) - a k (t)]. (2) 

The RF SQUID can be positioned near one of 
the antinodes of the vacuum current [e.g., see 
Figs. [21 (a) and (b) for k Q = 1] and can be biased to yield 
maximum coupling for any two of its eigenstates \i) and 
\j). The interaction matrix element between these lev- 
els represents their coupling strength with mode k and 
it is used to define the vacuum Rabi frequency gij — 
— (M cs / L S )I® k (z)(i\<&\j) /h. Moreover, when operating 
the system in the dispersive regime, the corresponding 
RF SQUID decay rates 7^ become 7?? = 7ij(#ij/(%) 2 , 
where Sij is the detuning between mode k and the transi- 
tion under consideration. The coupling #he, the effective 
decay rate 7^, and other relevant quantities have been 
calculated for both RF SQUIDs and PC qubits and typ- 
ical results are reported in Table I. 

A main application of the system illustrated above is 
the generation of single photons at frequency fk in a man- 
ner similar to a quantum-optical Raman scheme 0, ^| . 
After preparing the RF SQUID in level |g), the transi- 
tion |g) <-> |h) is driven by a classical excitation with Rabi 
frequency ^ g h ~ ghe and detuned by the amount S. The 
same transition is detuned from the resonator mode k by 
an amount A^> 5, resulting in a comparatively negligible 
coupling. On the other hand, the |h) <-> |e) transition is 
the only one coupled to mode k and it is also detuned by 
5 [see Figs. ^ (a) and (b)]. Choosing S ^> max[(] g h,^he], 
level |h) can be adiabatically eliminated thus lead- 

ing to the effective second-order Hamiltonian 

ffeff = S^|g>(g| + ^|e>(e|ata fc + 

+ % eff (|g)(e|a fc + |e)(g|a£) , (3) 

where g e ^ = (^ g h/^)#he is the effective Raman coupling. 
The first two terms at the r.h.s. of Eq. JSJ) are AC 
Zeeman shifts, while the last term describes an effec- 
tive anti-JC dynamics, inducing transitions within the 
{|g)|n), |e)|n + 1)} subspaces. The AC Zeeman shifts 
associated with the transition of interest {|g) |0), |e) |1)}, 
can be compensated by retuning the classical driving fre- 
quency. When the strong-coupling regime is reached, 
g e ff ^ max [ft c ,7hf], an effective 7r-pulse realizes a com- 
plete transfer of population from state |g)|0) to state 
|e)|l). This process leads to the creation of a microwave 
Fock state |1) inside the resonator that will leak out in 
a time ~ 1/k c . In the case of weak-coupling, the photon 
leaks to the outer world as soon as it is generated inside 
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the cavity, thereby realizing a deterministic single-photon 
source. Tailoring the photon pulse shape would require 
a time-dependent classical driving f2 g h (t) |8|. 

If the initial qubit-cavity state is (cos#|g) + 
e 1 ^ sin 0|e))|O), then the above Raman 7r-pulse would map 
it onto the cavity field (|e)|0) is a dark state of the anti- 
JC dynamics). In this way, we would be able to produce 
an outgoing field state cos#|l) + e %< ^ sin#|0). 

The proposed Raman scheme allows for generating sin- 
gle photons without initialization of the qubit in the ex- 
cited state |e) (|g)|0) is not a dark state of the anti-JC 
dynamics). Raman pulses are fast compare to STIRAP 
(adiabatic) techniques 7] and exploit the well defined 
qubit-cavity coupling |g, intrinsic in our scheme. 

Measurement schemes based on classical homodyn- 
ing are insufficient to resolve nonclassical field states. 
On the other hand, microwave single-photon detectors in 
mesoscopic systems do not exist to our knowledge. Here, 
we propose an MQHM technique as a means to measure 
relevant observables of weak quantum signals, even at the 
level of single photons. This can be implemented in three 
main steps. First, a signal (S) and a local oscillator (LO), 
characterized by the same angular frequency cos = ^lo, 
are coherently superposed at a suitably designed MBS 
cooled to base temperature [Fig. (a)]. Second, the mi- 
crowave fields at the MBS output ports are amplified 
at low temperatures by means of linear phase-insensitive 
amplifiers. Third, the amplified signals are then down- 
converted to DC currents / s?2 /4? proportional to the ener- 
gies of the input signals, via square-law detectors (and/or 
mixers) at room temperature 19]. The DC currents can 
be measured as voltages with an oscilloscope. Adequate 
manipulation of these measurements will lead to signal 
information with minimal noise background. 

The MBS is realized using a suitable four-port device: 
the hybrid ring [see Figs. El (a) and (b)]. The advanta- 
geous coplanar design proposed here can easily be scaled 
and preferably fabricated with low resistivity conductors. 
We now extend the classical theory of hybrid rings in 
Ref. yjj to the quantum regime by analogy with an op- 
tical beam splitter. With only the vacuum incident at 
ports two and four, and up to a global phase common 
to both input beams, the reduced quantum input-output 



TABLE I: Typical parameters and relevant quantities for RF 
SQUIDs and PC qubits, coupled to a 50 Q, resonator |22| . 
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FIG. 2: (Color online) (a) Sketch of the entire generation and 
measurement network (Cm, Ci, and C2' capacitive II network 
representing the resonator input port acting as a mirror; C p : 
parasitic capacitor at the cavity open-circuit ending; HR: hy- 
brid ring; PS: phase shifter; AT: attenuator), (b) Asymmet- 
ric CWG resonator with integrated HR. (c) On-chip antenna 
providing the classical driving Q g h. (d) Resonator-CWG cou- 
pling region, (e) The waves traveling around the HR interfere 
resulting in the plotted reflection and transmission amplitude 
patterns, (f) Isolation between ports three and one of the HR. 



relations of a lossless MBS are 
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Here, r and t are the complex, frequency-dependent re- 
flection and transmission coefficients, a s and the 
signal and LO port operators respectively. The latter 
is chosen to be a classical coherent field that is charac- 
terized by its complex amplitude a LO = |a LO | exp(i# r ), 
where \a LO \ is the norm of the field amplitude and 6 Y is 
its relative phase with respect to S. The numerical sim- 
ulations plotted in Fig. [21 (e) show that the MBS can be 
balanced over a broad bandwidth around the desired op- 
eration frequency i.e., r = t = l/y/2 (—3 dB), with 

outputs a 2 — (flg + ^lo)/ ^ an d a 4 = ( — a s + ^lo)/ 

At this point, standard optical quantum homodyning 
would require the use of photodetectors at each output 
port of the MBS. This would allow different measure- 
ments of photocurrents ^ p ,2/4? proportional to different 
realizations N of the observable ft — a) a. By comput- 
ing the difference — ^p,4 oc N2 — N4 = 2\a^o\Xo r , 
where the LO operators were replaced by their com- 
plex amplitudes, realizations of the quadrature Xq t = 
(age 2(9r + a s e~* 6>r )/2, enhanced by a factor 2|q;lo|, could 
be obtained. With the complete histogram of these mea- 
surements, all moment averages (Xg ), \/p G N, and 
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full reconstruction of the associated Wigner function via 
quantum tomography might be evaluated 5]. In ab- 
sence of microwave photodetectors, the MBS output sig- 
nals must go through linear amplifiers, square-law detec- 
tors and/or mixers before being measured at the oscillo- 
scope. These conditions impose severe restrictions in the 
measurement process with no counterpart in the optical 
regime and require additional theoretical considerations. 

The output signals, referred to the input, of linear 
phase-insensitive amplifiers can be written as Xa 2 — 
X a2 + X^ 2 and X~ a4 X a4 + X^ 4 , where X a = (a + a t )/2 
and X^ = (£ + i^)/2 The added noises at each 

arm, £ 2 /4> are random, uncorrelated, and characterized 
by (almost) the same noise temperature T n , leading to 
a mean photon number (n^) = = k^,T n /hwk- The 

difference between the measured currents produced by 
the square-law detectors, I s2 — ^s,4? is proportional to 
different realizations of the following observable 

it £ _ *t A , -t - i ;t; etc 
a 2 a 2 — a^a^ — a s a LO -+- & LO a s ' <?2S2 — S4S4 

, gt (gs+Ao) i g (4 + a Lo) 

^ t (-a s + a LO ) - (~4 + a Lo) /kn 

Repeating the measurement procedure, we can average 

(a 2 a 2 - a\a A ) = 2|a LO |(X^ r ), (6) 

given that (£2) = (£4) = and the reasonable assump- 
tion (n£ 2 ) ~ (% 4 )- Equation ® shows that the pro- 
posed MQHM allows the measurement of the enhanced 
mean value of the quadrature Xq t with negligible noise 
disturbance. This important physical quantity is sensi- 
tive to coherence: it is zero for any Fock state \n) and 
cos# r /2 for the superposition (|0) + \1))/V2. However, 
the method illustrated above does not lead to a mea- 
surement of (Xq ) since the amplifier noise shadows the 
information contained in the signal. Instead, we now 
propose to send both MBS outputs through mixers, us- 
ing the same calibrated local oscillator, and then evaluate 
the ensemble average of the measured products 

Here, we used (X^ 2 X^ 4 ) = (X^ 2 )(X^ 4 ) = and assumed 
(X^Q Gr ) to be known (after performing an adequate net- 
work calibration) or removable (e.g., via a modulation 
technique). Equation Q shows a remarkably simple way 
of measuring (Xq ) with minimal noise disturbance. 

A precise measurement of (Xq y ) and (Xq ) , as shown in 
Eqs. (0 and Q, requires random and uncorrelated noise, 
a sufficient number of repetition measurements, and ad- 
equate calibration, if the difference of (n^) ((%) ~ 21 for 



T n ~ 10 K) for the two amplifiers is not sufficiently small. 
It is known that the knowledge of (Xq t ) and (Xq ) pro- 
vides complete information about Gaussian states and a 
simple criterion for discriminating Fock states. Further- 
more, we conjecture here on the possibility of measuring 
(XI), Vp G N, V6> r e [0,2tt) [see PS in Fig. (a)] un- 
der the conditions described above, allowing a complete 
reconstruction of the Wigner function in the microwave 
domain. 

In conclusion, we proposed a new scheme for the de- 
terministic generation of intracavity and propagating mi- 
crowave Fock states or linear superpositions of them. We 
showed also how to realize MQHM for measuring first and 
second-order field quadrature moments. These proposals 
are essential tools for the implementation of quantum- 
optical CQED and linear optics in the microwave domain 
with superconducting devices on a chip. 

The authors thank C. M. Caves, S. M. Girvin, and 
D. C. Glattli for useful discussions. MM and WDO 
would like to acknowledge K. R. Brown, D. E. Oates, 
and Y. Nakamura for fruitful discussions and T. P. Or- 
lando and M. Gouker for their support. This work has 
been partially supported by the DFG through SFB 631. 
ES acknowledges EU support through the RESQ project. 



* Electronic address: Matteo.Mariantoni@wmi.badw.de 
[1] J.M. Raimond, M. Brune, and S. Haroche, Rev. Mod. 

Phys. 73, 565 (2001). 
[2] H. Mabuchi and A.C. Doherty, Science 298, 1372 (2002). 
[3] M.A. Nielsen and I.L. Chuang, Quantum Computation 
and Quantum Information, (Cambridge University Press, 
Cambridge, 2000). 
[4] D. Leibfried, R. Blatt, C. Monroe, and D. Wineland, Rev. 

Mod. Phys. 75, 281 (2003). 
[5] U. Leonhardt, Measuring the Quantum State of Light 

(Cambridge University Press, Cambridge, 1997). 
[6] A.I. Lvovsky et a/., Phys. Rev. Lett. 87, 050402 (2001). 
[7] A. Kuhn, M. Hennrich, and G. Rempe, Phys. Rev. Lett. 

89, 067901 (2002). 
[8] M. Keller et a/., Nature 431, 1075 (2004). 
[9] A. Badolato et al, Science 308, 1158 (2005). 
[10] A. Wallraff et a/., Nature 431, 162 (2004). 
[11] I. Chiorescu et a/., Nature 431, 159 (2004). 
[12] R. Movshovich et al, Phys. Rev. Lett. 65, 1419 (1990). 
[13] J. Gabelli et a/., Phys. Rev. Lett. 93, 056801 (2004). 
[14] M. Franca Santos, E. Solano, and R.L. de Matos Filho, 

Phys. Rev. Lett. 87, 093601 (2001). 
[15] C.-P. Yang, S.-I Chu, and S. Han, Phys. Rev. Lett. 92, 
117902 (2004). 

[16] K.V.R.M. Murali et a/., Phys. Rev. Lett. 93, 087003 
(2004). 

[17] Yu-Xi Liu, L.F. Wei, and F. Nori, Europhys. Lett. 67 (6), 
941 (2004). 

[18] J. Siewert et a/., e-print cond-mat/0509735 

[19] R.E. Collin, Foundations for Microwave Engineering, 2nd 

ed. (Wiley-IEEE Press, New Jersey, 2000). 
[20] M. Sarovar, Hsi-S. Goan, T.P Spiller, and G.J. Milburn, 



5 



Phys. Rev. A 72, 062327, (2005). and I « 435 nHm -1 respectively, and length D « 9.2 mm. 

[21] A. Blais et a/., Phys. Rev. A 69, 062320 (2004). [23] CM. Caves, Phys. Rev. D 26, 1817 (1982). 

[22] E.g., conductor width W = 118 /mi, slot 5 = 80 /xm, ca- 
pacitance and inductance per unit length c w 0.17 nFm -1 



